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Abstract 

We examine how the signature of the strange dibaryon resonances shows up in scattering am- 
phtudes and observables of the three-body KNN-ttYN (Y = S, A) system on the physical real 
energy axis. The so-called point method is applied to handle logarithmic singularities that ap- 
pear in solving the Alt-Grassberger-Sandhas equations for the real scattering energies. By taking 
two different kinds of models for the two-body KN-ttT, subsystem, both of which reproduce the 
available data equally well but give quite a different resonance-pole structure for A(1405), we also 
investigate whether the strange-dibaryon production reactions can be used for disentangling the 
nature of A(1405). 
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I. INTRODUCTION 



In recent years, the strange dibaryons have been studied actively as the simplest kaonic 
nuclei [l| in the three-body KNN-irYN system. A number of theoretical studies to search 
for the strange dibaryons have been performed with the variational method IJS and the 
Alt-Grassberger-Sandhas (AGS) equations IsMsl, employing the phenomenological potentials 
I2, ^jls] or the effective chiral Lagrangian [sljlsl,]?-^! for the meson-baryon and baryon-baryon 
interactions. All the studies support the existence of the strange dibaryons as resonance 
states in the energy region between the KNN and vrSA^ thresholds. However, the resonance 
energies predicted in those studies arestUl highly model-dependent. For example, the models 

give the resonance energies lower than those 



with energy independent potentials 

m 
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with energy dependent potentials 

In parallel with the theoretical works mentioned above, experimental searches for the 



lOj, the OBELIX 



strange dibaryons have also been done by the FINUDA Collaboration 
Collaboration 11 1, and the DISTO Collaboration [l^. Further data will become available 



_Collaboration [l^), and new 



16| experiments) and DA$NE 



from SPring-8 (LEPS Collaboration [13|) and GSI (FOP 
experiments are planned at J-PARC (E15 |l^ and E27 
(AMADEUS Collaboration ItI). 

In our previous works 7|, |8|], we have investigated a possible existence of the strange 
dibaryon resonances in the three-body KNN-tiTjN system. This has been achieved by 
searching for resonance poles of the three-body amplitudes in the complex energy plane, 
where the amplitudes are obtained by solving the coupled-channel AGS equations. There, 
two models, the energy-independent (E-indep) and the energy-dependent (E-dep) models, 
have been employed for the s-wave meson-baryon interactions, both of which are derived 
from the leading-order term of the effective chiral Lagrangian but those have different off- 
shell behavior. As a result, we have found one resonance pole of the strange dibaryon for the 
E-indep model while two for the E-dep model, which is summarized in Table [H This result 
indicates that off-shell behavior of the meson-baryon interactions of the two-body KN-nY 
subsystem is crucial for the resulting pole positions of the strange dibaryon resonances. 

Most of the theoretical studies have presented just pole positions of the strange dibaryon 
resonances. However, those are not a quantity that can be directly measured in experiments. 
In order to examine the existence of the strange dibaryons in connection with experiments. 
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(b) 



FIG. 1: Examples of the typical (a)kaon- and (b)photon-induced strange dibaryon production 
reactions. The strange dibaryon resonances would be produced in thick shaded boxes. 



TABLE I: Pole masses Mn of the strange dibaryon resonances obtained in our previous works 
See the text for the explanation on the E-indep and E-dep models. 



Re(Mij) (MeV) 



-Im(M/j) (MeV) 



E-indep model 
E-dep model 



2312-2326 
2354-2361 
2281-2303 



17-20 
17-23 
122-160 



one has to compute the cross sections of strange-dibaryon production reactions consistently 
in the same framework. The strange dibaryon resonances can be produced via, for example, 
kaon- and photon-induced reactions on light nuclei such as ^He and deuteron (Fig. [1]). Then 
the signal of the resonances would be observed in the invariant-mass and/or missing- mass 
distributions of the decay products. A couple of such studies have been performed by Koike- 
Harada and Yamagata et al. 19 1 on the basis of the optical potential approach. 

In this work, we examine how the signature of the strange dibaryons shows up in the ob- 
servables of the three-body reactions by applying our approach based on the coupled-channel 
AGS equations developed in Refs. 0, S]. It is well known that logarithmic singularities ap- 
pear when one solves the AGS equations for the break-up reactions at real scattering ener- 
gies. We handle those sin gula rities numerically by making use of the so-called point-method 



proposed by Schlessinger [22| and developed by Kamada et al. [5 



. With this method, we 



examine the behavior of the quasi-two-body amplitudes (the thick shaded boxes in Fig. [T]) of 
the KNN-ttYN system at real scattering energies between the KNN and tiTjN thresholds. 
Also, as a first step toward developing a model to compute reaction cross sections measured 
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at the facilities such as J-PARC and SPring-8 (e.g., the reactions in Fig. [T]), we examine the 
"transition probabihty" of a strange-dibaryon production reaction, {Yk)i=o+N — t- tt+S+A^, 
where {Yx)i=o is an "isobar" of KN states with isospin 1 = 0. 

In Sec. ini we explain the AGS equations for the three-body KNN-ttYN{Y = S, A) 
system and present the transition probability formula for break-up reactions. Then, we 
present the two-body meson-baryon interactions used in this work in Sec. IIIII The computed 
quasi-two-body amplitudes as well as transition probabilities for {Yk)i=o + N ^ tt + 1] + N 
are presented in Sec. IIVI Summary is given in Sec. |Vl A brief description of the point method 
is presented in Appendix. 



II. THREE-BODY EQUATIONS 

A. Alt-Grassberger-Sandhas equations 

Throughout this paper, we assume that the three-body processes take place via separable 
two-body interactions, which have the following form in the two-body center-of-mass (CM) 
frame, 

where g(^a)i {<li) is the cutoff factor of the two-body channel a{= jk) with a relative momentum 
Qi and isospin J; E is the total energy of the two-body system. In the three-body system, 
we define the two-body energy E as E = W — Ei{pi) with three-body energy W and the 
spectator particle energy Ei{pi), where pi is the relative momentum of the spectator particle 
i. The explicit forms of each two-body interaction are presented in detail in Sec. IIIII 

The assumption above implies that two-body subsystems in the three-body processes 
form an "isobar" and thus the processes can be described as a quasi-two-body scattering of 
the isobar and the spectator particle. The quasi-two-body amplitudes, X(^a)ii,ii3)iijiPi^'Pj'^ 



are then obtained by solving the AGS equations 20|, l2l[ | 



+ X] X^X] / dPnZi»)li,Mi"niPhPn,W) 

/ ^ \ / s:\ Til J 



X r(^),„n,(5),,m (W^ - En{Pn),Pn) X^S)j„n,{l3)pj{Pn,Pj, W) . (2) 
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TABLE II: Indices specifying the "isobars". 



Isobar 


Allowed isospin(s) 


Spectator particle 


three-body Fock space 


(Yk) = {KN2)AKNi) 

\ 1\ / \ -L/ 


0, 1 


Ni,N2 


\N1N2K) 


(y^) = (vrS) 

\ / \ / 


0, 1 


N 


liVSTr) ,|EiV7r) 

1 / ' 1 / 


(y^) = (vrA) 


1 


N 


liYAvr) ,|AiV7r) 

1 / ' 1 / 


{d) = (NN) 


1 


K 


\N1N2K) 


(N*) = (ttN) 


1/2, 3/2 


S 


\T.Ntt) ,\NT.tt) 


(N*) = (ttN) 


1/2 


A 


lAiVvr) ,|AfA7r) 


(dy) = m) 


1/2, 3/2 


TT 


\T.Ntt) ,\NT.tt) 


(dy) = (AN) 


1/2 


TT 


\KNtt) ,\NkT:) 



Here, (a)/ denotes the isobar formed by a two-particle pair a with isospin J; the subscripts 
n represent the spectator particles. The notations for the isobars are summarized in 
Table HI 

The driving term Z(^a)ii,(i3)j,jiPi,Pj',W) describes a particle-exchange potential given by 
[see Fig. |2] (a) for the kinematics] 

W) = + (3) 

where Ei{pi) and Ej{pj) are the energies of the spectator particles i and j, respec- 
tively; Ek{pk) with pk = —pi — pj is the energy of the exchange particle k; qi {(jj) is 
the relative momentum between the exchange-particle and the spectator-particle j (i). 
In the nonrelativistic kinematics, we have En{pn) = f^n + pi/C^i^n) {n = i,j,k) and 
qi,j = {mk,iPj,k - mj^kPk,i)/{mj,k + nik^i). The s-wave projection of Z(^a)ii,{p)pj{PuPj]W) 
is given by 

1 

Zi»)ii,WidiPi^Pj'^^) = 2 J ^(cos^)^(«)/»,(/3),/i(Pi)Pi; (4) 
with cos 9 = Pi - Pj. 

The isobar propagator, T<^a)ii,{p)ii (W ~ Ei{pi),pi) as illustrated in Fig. [2] (b), is given in 
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(a) 



(b) 



FIG. 2: (a) One particle exchange interaction Z(^^-^ji(^p-^^^j{pi,pj^W). (b) Isobar propagator 
the nonrelativistic kinematics by solving the following Lippmann-Schwinger equations, 

T{a)ii,{l3)ii{^ - Ei{pi),Pi) = \{a)ii,{l3)ii 

+|Jy ''^'Hv - Em - £,.(p..*)'™'-'^"''''^ " ^^^'"'^ ■ 

(5) 

Here, Ej^^pi, qj) is the energy of the interacting pair (jk), Ejk{pi, qi) = rrij + rrik+pf /2{mj + 
mfc) + gf/2/ij with the reduced mass defined as /ij = mjimk/ {mj + m^). 

After taking antisymmetrization for the two-nucleon states in the three-body processes, 
the AGS equations ([2]) are formally written as (suppressing all indices other than that of 
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the isobars) 



(x 



{Yk),{Yk) 



X 



{Y.)XYk) 



X. 



WAYk) 



X 



{N*),{Yk) 



( ^{YkUYk)^ 



\ 



<d),{YK) 






^Z(YK),(YK)T(YK),iYK) Z(YK),iYK)HYK),iY^) Z (^Yk) ,id)^d) 



Z{dUYK)r(YK),{YK) Z(d)XYK)nYK)XY.) 

Z{N*),(Y^)r(Y^),(YK) Z(^N*),(Y^)T(Y^),{Y^) 

\ Zi^dy),{Y^)r(Y^),{YK) Z(_dy),(Y^)T{Y^),{Y^) 



X 



X 



{Y.),{Yk) 



X, 



(d),iYR 



X 



(N*),{Yk) 

\X(dy),{YK)/ 



id) \ 



Z^N*),(dy)ndy),{dy) 

Z(dy),{N*)T(N*),{N'') / 



(6) 



B. Break-up reactions 

In this subsection, we present formulas for computing transition probabihty of the quasi- 
2-body to 3-body reaction, {Yk)i=o + A^— t-tt + S + A^. For this purpose, we first need to 
define the amphtudes of the {Yk)i=o + A^— t-tt + S + A^ reaction. This is because within 
our formulation the well-defined amplitudes are of the 3-body to 3-body scatterings, where 
all the external particles are stable against strong interactions. The relevant amplitude here 
is of the {K + N) + N ^ (^i^)7=o + iV— t-tt + S + A^ reaction, which is given in a concise 
notation as 

Tjrj:N^(KN)N = X] ^^9la)iT{a)ii,{'y)ii^{-f)ii,{YK)i=oNT{YK)i=oN,{YK)i=oNg{YK)i=o^ (^) 
{a)i=-KT,N (7) / 

where the summation of (a)z is taken for all possible combinations of ttJ^N. Now let us 
consider the isobar (Yx)i=o as an actual resonance state of the two-body reactions. (Note 
that we originally introduced notion of the isobars just for the sake of convenience in our 
formulation and did not take them as actual resonances.) Near a resonance pole of the 



isobar propagator T"(y^)j^o,(y^)^^o(-E'), the two-body amplitude for KNi=q — )■ KNj=o can be 
approximated as 

(E) 

^ 9(Yj,)j^o E-M + iT/2 ^(^^)^=o 
= 9*kn^Yk e - M + iT /2^^^^^^' 

where M — iV/2 is the resonance pole position of T(y^)j^^o.(>K)/=o(-^) ^^"^ ^{Yk)i=o the 
residue of T(j^)i^o,{Yk)i=o at the pole. Also, gYK*-KN = \/ R{Yk)i=o9{Yk)i=o [Qkn^Yk = 
\/ R{Yk)i=o9{Yk)i-o\ '^^^ interpreted as a vertex function for the process KNj=o — (Yk)i=o 
[{Yk)i=o KNi=q]. Within this approximation, the three-body amplitude can be written 
as 

{a)i=TTT.N (7) / 

X ^(yK),=„^,(ix)z=o^(^ - ^7v(p7v),P7v)c/(yK),=o 

(a)i=7rSAr (7) / 



X 



R{Yk),^oG(j^)^^^{W - EN{pN)jN)gYK^KN. (9) 



where G(Yk)j-(,(W ~ En{pn),'Pn) is the (Yk)i=o resonance propagator in the existence of a 
spectator nucleon with momentum p^. From Eq. (|9]), it is reasonable to define T-matrix of 
{Yk)i=o + N^7r + J: + Nas 

Tt,t.N'^(Yk)i=oN = y^,y^,9*a)jMii,{'-t)ii^i'y)ii,{yK)i=oN\/RiYK)i=o- (10) 

{a)i=TTT,N (7) / 

The s-wave projection of the scattering amplitudes for the {Yk)i=o + N^7i + J] + N 
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reaction are then given by 

I 

X {| [[tt ® S](y^), ® iV]r) g{Y^)j{qN)T(Y^)jN,{YK)iN {W - En{pn)) X^Yk)iN,{Yk),=on{.Pn .p'n. W) 

+ \[[k® S](y^), ® N]y) g{Y^)MN)r{Y^)jN,{Y^)iN [W - En{Pn)) X(y^)jN,(Yk)i^MPN ,p'n ^ W) 

+ |[[7r ® N\n*)i ® S]r) g(N*)M^)nN')iJi,{N*)iJi {W - E^{py^)) X^^n')iK,{Yk)i=on{p^^p'n^ ^) 
+ I [[S N\dy)i ® vr]r) g{dy)A(l^)ndy)i^,{dy),^ {W - E^{p^)) X^dy)jn,iYK)i=oN{Pn,p'N: W)} 

x([(F;,),=o®ivviv/iwi;:, (11) 

where \ [A B]a ® C]b) with (ABC) = (ttSA^) and | [{Yk)i=o ® A^]^) are the spin-isospin wave 
functions of the final and initial states, and Xi^a)ii,{i3)^ij{PiP'i W) is the s-wave projection of 
the quasi-two-body amplitudes given in Eqs. (Ej) and ([6]). The momentum gs, ^tt, Pt. and p^,- 
are function of qn and pW, i.e., qT.{qN,PN), qn{qN,PN), PT.{qN,PN) and Pn{qN,PN)- 

Using Eq. ( ITT]) , we define the transition probability of {Yk)i=o + N-^7r + T, + Na.s 
follows, 

(12) 

III. MODEL OF TWO-BODY INTERACTIONS 

Now we present explicit forms of the two-body interactions [Eq. ([T])] used in this work. 
We first consider the meson-baryon interactions (Sec. IIII Al) . and then consider the baryon- 
baryon interactions (Sec. IIII Bl) . In this section, we suppress indices of the spectator. 



A. Meson-baryon interaction 

As done in our earlier works j?, 8|, we consider two kinds of models for the s-wave meson- 
baryon interactions, which are called the E-indep and E-dep models, respectively. The 
explicit forms are given by 

^E-indcp / / 
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(13) 



for the E-indep model, and by 

= 32^2^2 ^/^^^'' •^("^^^^'^^(^^^^^^' ^^^^ 

for the E-dep model. Here, rria (M^) is the meson (baryon) mass of the channel a; q' 
(g) is the magnitude of relative momentum of the channel a in the two-body CM 
frame; is the pion decay constant; the coupling coefficients C(^a)i(i3)i are summarized in 
Table UTTl As for the cutoff factors g(a)i{Q'), "we employ the dipole form with the cutoff A(^a)n 



TABLE III: The coupling coefficients C'(a)j(/3)j- Note that Ci^a)i{p)i = 



(a,/?) 


Total Isospin / 




{KN, KN) 





6 


{KN, vrS) 





-Vd 


(vrS, vrS) 





8 


{KN, KN) 


1 


2 


{KN, vrS) 


1 


-2 


{KN, vrA) 


1 


-V6 


(vrS, vrS) 


1 


4 


(7rS,7rA) 


1 





(vrA, vrA) 


1 





(vr Af, vriV) 


1/2 


4 


{irN, ttN) 


3/2 


-2 



It is noted that except for the cutoff factors, both of the above potentials [Eqs. ([I3 



and ([H])] are derived from the so-called Weinberg- Tomozawa term 2J,|25| that is the leading- 
order term of the effective chiral Lagrangian, 



WT = TT^tr(^B7^[[</., ^Ab]), (15) 



8F2 

with ipB (0) being the octet baryon (pseudoscalar meson) field. From this Lagrangian, the 
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s-wave potential is given by 



t^WT 



327r2F2 ^uj^{q')ujfs{q) V 2E„ iq')2E(s (g) 

X [oo^iq') + - M„ + uj^{q) + ^^^(g) - M^] , (16) 

where 0Ja{q') [Ea{q')] is the meson [baryon] energy of the channel a. We then obtain the 
E-indep potential f lT3|) from Eq. f|T6|) by assuming |g '| ^ rria, Ma and |g| ^ mp.Mp. On 
the other hand, the E-dep potential ( fT4l) are given by first replacing Ua{q') + Ea{q') and 
oopiq) + Ep{q) in the brackets of Eq. f ITB]) with the on-shell two-body scattering energy i?, 
which is now considered to be an independent variable, and then assuming \q'\ <^ ma,Ma 
and |g| ^ m^, M^. The replacement with the on-shell two-body scattering energy in deriving 



26| 



the E-dep potential corresponds to the so-called "on-shell factorization'' 

As already seen in Sec. [Tll we take the non-relativistic kinematics for the numerical 
calculations. This is because of a problem inherent in the use of energy- dependent two-body 
potentials for the three-body calculations with the relativistic kinematics. If the relativistic 
kinematics are used, the total energy of the two-body subsystem can become pure imaginary 
fo.. ,a,,e spectator ™o.e„ta fl. Howeve. sucK a d,mcu.ty does not appear as far as one 
uses the non-relativistic kinematics. 

Parameters of the two-body potentials are the cutoffs ^{a)i- We determine the cutoffs by 
fitting the / = vrS invariant mass distributions of the K~p — )■ ttttttS reaction and the KN 
reaction cross sections. Results of the fit for the E-indep and E-dep models are presented 
in Fig. |3] and Fig. HI respectively. There, the results are shown as bands because we have 
determined the cutoffs only up to certain ranges within which the computed cross sections 
are consistent with the experimental errors. The fitted values of the cutoffs are listed in 
Table IVl 

TABLE IV: Cutoff parameters of KN-nY interaction. 





A{y^),=o (MeV) A(y^),^„ (MeV) A(y,,),^, (MeV) A(y^ 


=7rS)/^l 


(MeV) AiY^=^A)r^, (MeV) 


E-indep 


975-1000 675-725 920 


960 


640 


E-dep 


975-1000 675-725 725 


725 


725 



In Fig. [5l we present the resonance pole positions of the KN s-wave scattering amplitudes 
in the complex energy plane between the KN and ttE threshold energies. We find that 
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(a) 



\\ 






100 150 200 250 50 100 150 200 




(e) 



200 250 50 100 150 200 250 



FIG. 3: Results of the fit with the E-indep model, (a) I = vrS invariant mass distributions of 
K^p — vrvrvrS; total cross sections of (b) K~p — )■ K~p, (c) K~p — )• tt'^S^, (d) K~p — >• 7r~S+, and 



{e) K p — 7- vr S . Data are from Refs. 
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FIG. 4: Results of the fit with the E-dep model, (a) / = vrS invariant mass distributions of 
K^p — 7- TrvrvrS; total cross sections of (b) K~p — t- K~p, (c) — t- 7r~^S^, (d) K~p — )■ 7r~5]~^, and 



(e) K p — 7- vr S . Data are from Refs. 



27 



^1 



the E-indep model has a single pole corresponding to A(1405) in the KN physical and vrE 
unphysical sheet [FiglS](a)], while the E-dep model has two poles in the same sheet [FigO^b)]. 
The analytic structure of the amplitudes in the E-dep model is similar to that obtained with 



the chiral unitary model 



33| 
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1350 1400 
Re[E] (MeV) 




1350 1400 
Re[E] (MeV) 



FIG. 5: The S = —1, = 1/2 KN s-wave amplitude on complex energy plane in (a) the E- 
indep model and (b) the E-dep model. The cutoff parameters are (A^y^)^^^ , A^y^)^^^, ) = (1000, 700) 
MeV. 



As for the cutoffs 
scattering lengths 



with a = TiN, we have determined them by fitting the Su and 5*31 vrA^ 
34| . The resulting values are summarized in Table |Vl 



TABLE V: Cutoffs of the vrA^ interactions. 





A(7V.),.,/,(MeV) 


A(7V.),.3/.(MeV) 


E-indep model 


400 


400 


E-dep model 


400 


400 



B. Baryon-baryon interactions 

As for the s-wave A^A^ interactions, we take the following form [3], 

q) = A-KCRgR{q')gR{q) - A-KCAgA{q')gA{q)- (17) 

Here, Cr {Ca) is the coupling strength of the repulsive (attractive) potential; The form 
factors gR,A{q) are defined by gR,A{q) = ^R,A^/{q^ + ^r,a^) with Is^r^a being the cutoff 
parameters of the A^A^ interactions. The coupling strengths Cr^a and the cutoff parameters 



A^^^ are determined by fitting the ^Sq phase shifts |35| (see Fig. [6] for the result of the fit). 



The resulting values of the parameters are summarized in Table 



As for the s-wave YN interactions, we follow the form given in Ref. 36 1 



C 

Vic^hmii^'^q) = -4vr^|^^(/ia/i/3A(„),A(^)J"^/^^(„),(g')m,(g)- (18) 
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TABLE VI: Parameters of NN interaction. 



Ai?(MeV) 


AA(MeV) 


C/j(MeV fm^) 


CA(MeV fm^) 


1215 


352 


5.05 


5.84 




50 100 150 200 250 300 



E|ab (MeV) 



FIG. 6: Phase shifts of A^A^ scattering for the ^Sq state. The soUd hne shows the phase shift with 
our model, and triangles show the phase shifts with the model of Ref. j35l ]. 



Here, C(^a)i{i3)i are the coupling constants given by C(^zn)i{sn)i = 0.83, C(^ej^)j(^an)i = 0.56, 
and C(^AN)i{AN)i = 0.49 for / = 1/2, and C^y:n)i(en)i = —0.29 for / = 3/2; /i^ is the reduced 
mass for the YN system; the form factor g{a)i{<l) is defined by g[a)j{(l) = ^fa)i/(l'^ + ^W/-'' 
the cutoff parameters A(o)^ are given by A(sAr)^ = 251 MeV and A(AAr)^ = 262 MeV. 



IV. RESULTS AND DISCUSSION 



A. Quasi-two-body scatterings 

Now we present the partial-wave quasi-two-body amplitudes at the real scattering en- 
ergies W, X(^a)ii,{i3)j,j{Pi,Pj,W), which are obtained by solving the coupled-channel AGS 
equation (|6]) and using the point method explained in Appendix. In FigJTJ we present the 
absolute square of the amplitudes, \X(^YK)i=oN,(YK)i=oN{PiyPjy'^)\'^^ whose initial- and final- 
state isobars are (Yk) with the isospin / = 0. Here we plot the results of the E-indep (E-dep) 
model as solid (dashed) curves. Also, we plot the amplitudes with two different cases of the 
off-shell momentum for each model, with pi = pj = 150 MeV for thick curves and with 



14 




TiEN 2300 2350 KNN 

Wc.m. (MeV) 

FIG. 7: W dependence of \X(^YK)i=oN,{YK)i=oNiPi^Pjj^)\'^ ■ Each curve is: (solid curves) the E- 
indep model; (dashed curves) the E-dep model; (thick curves) pi = pj = 150 MeV; (thin curves) 
Pi = Pj = 100 MeV. The cutoff parameters are taken to be (A(y^)j-^y, A(y^)j-^g) = (1000, 700) MeV. 

Pi = Pj = 100 MeV for thin curves, to examine the momentum dependence of the ampli- 
tudes. We find both models have a bump between the KNN and nUN threshold energies: 
W ~ 2305 MeV for the E-indep model and W ~ 2340 MeV for the E-dep model, both of 
which are close to the resonance pole masses Mr with — Im(M/j) ~ 20 MeV (see Table [T]). 
Furthermore, the positions of the bumps are independent of the momentum, and thus we 
can conclude that these bumps are actually produced by the strange dibaryon resonances. 
On the other hand, in the E-dep model, another strange dibaryon with — Im(MR) ~ 100 
MeV barely affects the amplitude on the physical real energy axis. This is consistent with 
the fact that normally resonances with large widths cannot produce a sharp peak in the 
absolute square of the amplitudes or cross sections. In Fig. El we show the W dependence 
of the amplitudes with different final states. We observe that the bumps due to the strange 
dibaryon resonances appear at almost the same W regardless of the final quasi-two-body 
states, as it should be. 

Next we present the contributions of each one-particle-exchange mechanism Z to the am- 
plitude \X{YK)i^oN,{YK)i=oN{Pi,Pj, W^)p with Pi = Pj = 100 MeV (Fig.|9]). Here the solid curve 
in Fig. M^a) [Fig. IHt^b)] is same as the thin-solid (thin-dashed) curve in Fig. [3 If the baryon- 
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FIG. 8: W dependence of |X(c,)^i,yj^(/=o)Ar(PMPi> (a) (a)/ = {Yt,)i=o and i = N, (b) (a)/ = 

((i)/=i (repulsive) and i = K, (c) (a)/ = (A^*)/=i/2 and i = S, and (d) (a)/ = ((iy)/=i/2 and i = tt. 
The meaning of each curve and the cutoff parameters are taken to be the same as in Fig. [71 

exchange (meson-exchange) Z-potentials are switched off in the rescattering processes, then 
the sohd curves in Fig. [9] are turned into the dashed (dotted) curves. Contributions of the 
meson-exchange processes seem to be crucial for producing the similar bump structure to 
the full amplitudes, while those of the baryon-exchange processes do not. However, we also 
observe that rescattering effects including both the meson- and baryon-exchange processes, 
which are required by the three-body unitarity, amplify the magnitude of the scattering 
amplitudes significantly, indicating the importance of maintaining the three-body unitarity 
exactly in searching for the evidence of the strange dibaryon resonances. 



B. Transition probability for the break-up {Yk)i=o + N^tt + Tj + N reaction. 

Finally, we investigate the energy dependence of the transition probability, w{pn,W) 
defined in Eq. f|T2|) . for the (Yr-)/=o + A^— j-tt + E + A^ break-up reaction. In Fig. [TOl we 
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FIG. 9: Contributions of one-particle-exchange processes to \X(Yj^-j^^^j^^(^YK)i=oNiPi^Pjj^)\'^- The 
figures are for (a) the E-indep model and (b) the E-dep model. Each curve is: (solid curves) full 
results; (dashed curves) baryon-exchange processes only; (dotted curves) meson-exchange processes 
only. The momentum [cutoff parameters] are fixed as pi = pj = 100 MeV [(A(-y^^,)^^jj, A(y^)^^Q) = 
(1000,700) MeV]. 

present w{pn, W) for = 100 MeV and p^r = 150 MeV using the same values of parameters 
as used in Fig. [71 We again find that the position of the bumps in w{pn, W) are independent 
of the momentum p^ of the initial {Yk)N channel, implying that the bumps originate from 
the strange dibaryon resonances. The E-indep and E-dep models are found to produce 
quite different energy dependence on the transition probabilities; those differences would 
be large enough to be detected by experiments. Because this difference is closely related 
to the different nature of A(1405) between the two models as shown in Fig. the strange 
dibaryon production reactions would also provide critical information on the dynamical 
origin of A(1405). 

Next we examine the cutoff parameter dependence on the transition probability w{pn, W) 
(Fig. [11]). The bands are given by varying the values of A(^Yk)i=o ^Ci'7r)/=o within the 
allowed range listed in Table IIVI We see that the signal of the strange dibaryon resonances 
remains to be observed in the transition probability within the allowed range of Aiyi^-^j^q 
and A{yO/=o- 

Finally, we examine the contribution of each reaction process to the transition prob- 
ability (Fig. [T2|) . As can be seen in Eq. ( ITT]) , the reaction processes consist of the 
quasi two-body processes characterized by the amplitudes X(Yj^)jN^(^Yk)i=oN, ^{y^)iN,{Yk)i=oN-i 
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FIG. 10: Total transition probability w{pN, W) for {Yk)i=o + iV-;>7r + E + iV. The meaning of 
each curve and the cutoff parameters are taken to be same as in Fig. [71 
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FIG. 11: Cutoff dependence on the transition probability for the (Yk)i=o + N — > tt + T, + N 
reaction. Panel (a): the E-indep model; Panel (b) the E-dep model. The bands of transition 
probability are produced by varying values of A(^YK)i=a ^(Yn)i=o allowed range listed in 

Table ITVl The initial nucleon momentum is set to pjsf = 100 MeV. 

X(N*)iJ:,{yK)i=oN, and X^dy)i7r,{YK)i=oN- We find that the X^Yk)i=oN,{Yk)i=oN process has dom- 
inant contribution about 85% to the transition probabiUty, while others have rather small 
contributions: about 5% is from X(^dy)in,(YK)i=oN^ ^^^^ than 1% is from X(y^)jN,{Yk)i=oN 
and X(^n'-)iJ:,(Yk)i=oN- 
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FIG. 12: Contribution of each quasi-two-body process to the {Yk)i=o + A^— s- vr + E + iV transition 
probabihty. The figures are for (a) the E-indep model and (b) the E-dep model. Each curve 
is: (solid curve) Full results; (dashed curve) X(^Yk)iN,{Yk)i=oN process only; (dashed-two-dotted 
curve) X(^Yr:)iN,{YK)i=oN process only; (dotted curve) ^{Ar*)^E,{yA')/=oAf Process only; (dashed-dotted 
curve) ^{dy)iTT,{YK)i^oN process only. The cutoff parameters A(^Yk)i=o ^{Yn)i=o taken to be 
1000 MeV and 700 MeV, respectively, and the initial nucleon momentum is set to pat = 100 MeV. 

V. SUMMARY 

Within the framework of the coupled- channel AGS equations, we have examined how the 
signature of the strange dibaryon resonances in the three-body KNN-irYN system shows 
up in the scattering amplitudes and transition probabilities on the physical real energy axis. 
The logarithmic singularities appearing in solving the AGS equations for the real scattering 
energies have been successfully handled by making use of the point method. Two different 
kinds of models, the E-indep and E-dep models, have been considered for the two-body 
KN-TiT, subsystem to investigate whether the strange-dibaryon production reactions can be 
used for disentangling the nature of the two-body KN-ttTj system with A(1405). 

We have found that within our model, a clear bump produced by strange dibaryon 
resonances appear in the quasi-two-body scattering amplitudes ^(o)ji,{yK)/=o^(^) ^^"^ 
{Yk)i=o + N ^ tx + + N transition probabilities in the energy region between the KNN 
and vrSA^ thresholds, which strongly suggests that the clear signals of strange dibaryon reso- 
nances should be detected by measuring of vrSA^ invariant mass distributions at the relevant 
energies. We have also found that the E-indep and E-dep models produce quite different 
energy dependence on X(^a)ii,{YK)i=oNiW) and {Yk)i=o + A^— t-tt + S + A^ transition proba- 
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bilities; those differences would be large enough to be detected by experiments. Within our 
framework, this difference originates from the different nature of A (1405) between the two 
models as shown in Fig. El and thus the strange dibaryon production reactions would also 
be helpful to reveal the dynamical origin of A(1405). 

It is for the first time that the break-up (Yk)i=o + A^— t-tt + S + A^ transition proba- 
bilities are computed within the fully coupled-channel AGS equations. As a next step, we 
will further account for initial state interactions and develop a technique to make practical 
calculations of "actual" cross sections of kaon- and photon-induced strange dibaryon pro- 
duction reactions shown in Fig. [H which will be measured at experimental facilities such as 
J-PARC and SPring-8. This will be discussed elsewhere. 
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Appendix A: Brief description of the point method 



The s-wave projection of the particle-exchange potential, Z(^a)ii,ii3)i'jiPiyPjy'^) [Eq- 
contains the following logarithm. 



In 



n2 „2 



2mfc 



2mi 2mj Sm^ 



(Al) 



For real W with W > M, this logarithm becomes singular at momentum {pi,Pj) satisfying 



W -M 



. 



2mj 2mj ^rrik 
The singularities appear as a "moon-shape" in the prPj plane as illustrated in Fig. 



(A2) 
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FIG. 13: The moon-shape singularities. Sohd curve shows the momentum {pi,Pj) where 
^{a)ji,{i3)j,jiPijPjjW) has logarithmic singularity. 

As a practical technique to handle the moon-shape singularities in solving the scatte ring 
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equations ([6]), we have employed the point method, which is proposed by Schlessinger 
and developed by Kamada et al. 23|. We briefly explain the method in the following. 

The point method is an extrapolation technique of functions. With this technique, one 
can evaluate the value of a function X{W) of real W from X{W + iei), where (e = 1, 2, ...) 
is a series of positive finites that converges to zero, using the following formulas. 



X{W) = lim 



ai(^-^i) 
1+- 



1 



lim 

e-5.0 



X{W + iei) ai(e - ei) a2(e - ea) 



1+ 



1+ 



with 



1 



ai 



o-i-iiei+i — ez-i) a;_2(ei+i — e;_i) 



1- 



(A3) 



1+ 



Q - e/+i 

To illustrate how we get scattering amp 
the formulas above to the Amado model 



1+ 1-[X{W + tei)/X{W + iei+i)] 

itudes X(^a)ii,{^)j,j{Pi,Pj, W) for real W, we apply 



(A4) 



20[ |. a simple model for three-boson scatterings. 



The AGS equations for the s-wave scattering of a boson b and a two-6 bound-state d, bd — )■ bd, 
are given by 



X(p',po, W) = 2Z{p\p, W) + 2 p^dpZip',p, W)t{p, W)X{p,po, W). 



(A5) 
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In the bd CM system, the driving term Z{p' ,p,W) and the two-body propagator t{p,W) 
are expressed as 



Z{p',p,W) = ^ j'^dx-^ 



9o 90 



W'-f-f-^+^e ^^^^ 

2m 2m 2m 



T-\p, W) = [E^ip, W') + B + ie] 

9\k) 



{E2{p, W') + B + ie) j k' 



dk- 



(A7) 



(^ + 3 {Mp,W')-^ + ^e), 

Here, g{q) = 90/ {q"^ + f3^) is the form factor for d — t- bb, which is normahzed as 
/ k'^dkg'^{k)/{B + = 1; B is binding energy of d, and E2{p,W) = W - 'ip^/{Am) is 
two-body scattering energy. We solve these AGS equations by setting h = 2m = 1, B = 1.5, 
13 = 5, and W = 1. 

If one tries to solve Eq. flA5|) for a real W, the momentum integral path crosses the singu- 
larities of the Z potential and thus the resulting amplitude X{p',pQ, W) does not converge. 
On the other hand, one can have convergent solutions of Eq. flASP without any problems for 
complex energies W + iei with positive finites e^. Therefore, we first compute the amplitude 
X for several complex energies, and then make an extrapolation to X{W) using Eqs. ( ]A3p 
and ( ]A4I) . For practical computations, we use five e^'s, 

ei = 0.05x1 (/ = 1,2,...,5). (A8) 



In Fig. [HI we show the p dependence of X{p, pq, W) for = 1 and po = Am{W + -B)/3. 
The sohd (dashed) curve represents the real (imaginary) part of the amplitude X{p,pQ,W) 
extrapolated using the point method. In the same figure, we also present the amplitude 



obtained by the spline interpolation method 37| as a comparison. 



The scattering amplitude X{W) for the KNN-tiYN system studied in this work is ex- 
trapolated from the amplitude X{W + iei) at e/ = 10 x / (MeV) for / = 1, 2, . . . , 5, using 
Eqs. dAH) and f lAi]) . 
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